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ABSTRACT 


Thermals caused by convection in the lower atmosphere are commonly used by birds and 
glider pilots to extend flight duration, increase cross-country speed, improve range, or simply to 
conserve energy. Uninhabited Aerial Vehicles (UAVs) can also increase performance and reduce 
energy consumption by exploiting atmospheric convection. An autonomous soaring research project 
was conducted at the NASA Dryden Flight Research Center to evaluate the concept through flight 
test of an electric -powered motor-glider with a wingspan of 4.27 m (14 ft). The UAV’s commercial 
autopilot software was modified to include outer-loop soaring guidance and control. The aircraft 
total energy state was used to detect and soar within thermals. Estimated thermal size and position 
were used to calculate guidance commands for soaring flight. Results from a total of 23 thermal 
encounters show good performance of the guidance and control algorithms to autonomously detect 
and exploit thermals. The UAV had an average climb of 172 m (567 ft) during these encounters. 

NOMENCLATURE 

a equatorial radius of the earth, m 

dt queue sample time step, s 

H 

E total specific energy of the aircraft, m 

E specific energy rate, m/s 

E motor specific energy rate due to the motor, m/s 

E post flight estimated specific energy rate, m/s 
E specific energy acceleration, m/s 2 

E s scaled specific energy acceleration, m/s 2 

e p controller position error, m 

e„ Jim „ n mean sum squared velocity error, m 2 /s 2 

e v controller velocity error, m/s 

e w weighted velocity error, m/s 

e sumsq mean sum squared velocity error for perturbed thermal radius, m 2 /s 2 

e w weighted velocity error for perturbed thermal radius, m/s 

g gravitational constant of Earth, m/s 2 

H velocity influence vector 

h gst estimated aircraft altitude, m 

i queue index 

*z 

k current measurement index 

energy acceleration gain, deg*s/m 
k ? position error gain, deg/m* s 

velocity error gain, deg/m 
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learning rate, 1/s 2 

latitude of aircraft, rad 

latitude of X Y reference frame, rad 

longitude of aircraft, rad 

longitude of X Y reference frame, rad 

number of rows used in drift calculation 

number of rows in the queue 

aircraft mass, kg 

aircraft position, m 

estimated thermal position, m 

thermal position calculated from the first 20 entries in the queue, m 
thermal position calculated from the last 20 entries in the queue, m 
aircraft X position, m 
aircraft Y position, m 
static pressure, Pa 

queue of aircraft position and shifted energy rate 

queue of aircraft position and energy rate 

queue of drift corrected aircraft position and energy rate 

radius of commanded circular flight path, m 

estimated thermal radius, m 

perturbed thermal radius, m 

mode logic state 

distance between the aircraft and the thermal center, m 

Laplace transform variable 

thrust from the motor, N 

maximum thrust from the motor, N 

one-second delayed throttle command, fraction 

trim throttle needed for level flight, fraction 

true airspeed, m/s 

estimated thermal drift velocity, m/s 

sink rate of the air between thermals (environment sink rate), m/s 
estimated aircraft velocity rate caused by the motor, m/s 2 
predicted thermal vertical velocity, m/s 
estimated thermal vertical velocity, m/s 

predicted thermal vertical velocity for perturbed thermal radius, m/s 


dy thermal radius stepsize, m 

A perturbed thermal radius error gradient, m/s 2 

t \f cmd turn rate command, deg/s 

\j/ steady state turn rate, deg/s 

INTRODUCTION 

A novel way to improve the range, duration, or cross-country speed of an autonomous aircraft 
is to use buoyant plumes of air called thermals. Thermals occur when the air near the ground 
becomes less dense than the surrounding air as a result of heating or humidity changes at the 
Earth’s surface. Thermal height varies greatly depending on climate and season but can typically 
fall between 500 m (1600 ft) and 3700 m (12 000 ft) above ground level. Sailplane pilots primarily 
rely on thermals to travel great distances and soar for long periods of time (ref. 1). Large birds, such 
as hawks and buzzards, have been observed to soar extensively while searching for food and vary 
their migration path and flight times in response to thermal activity (refs. 2-4) and frigatebirds are 
known to soar continuously day and night using thermals over tropical ocean waters (ref. 5). 

Many unmanned aerial vehicles (UAVs) have similar sizes and wing loadings to soaring 
birds and manned sailplanes. Mission profiles that could allow small UAVs to take advantage 
of thermals include remote sensing, surveillance, atmospheric research, communications, force 
protection, forest-fire monitoring, land management, and border control. A study using a simple 
UAV simulation with thermals calculated from meteorological data from Desert Rock, Nevada, 
found that a UAV with a nominal endurance of 2 hours can gain up to 12 hours of flight time by 
using thermals (ref. 6). 

Many strategies for extracting energy from thermals and other atmospheric phenomena have 
been published for glider pilots (refs. 1 and 7). Autonomous soaring for UAVs is an emerging field 
of study. Boslough (ref. 8) used piloted flight tests with an instrumented model glider to investigate 
specific maneuvers flown in the presence of wind shear and proposed that flocks of soaring UAVs 
could radically increase the efficiency of powered flight. Patel showed through simulation that 
significant energy savings can be achieved for UAVs by performing pitch maneuvers in turbulent 
air (ref. 9). Autonomous thermal soaring UAVs were first proposed by Wharington (refs. 10 and 11) 
in 1998 as a method to extend UAV performance and possibly to enable perpetual flight. Optimal 
guidance algorithms were developed using reinforcement learning and a neural-based thermal 
locator to autonomously detect and utilize thermals. Results showed that both simple heuristics and 
reinforcement learning could be used along with the thermal locator to improve UAV performance, 
but the reinforcement learning algorithms were too computationally intensive for real-time use. 
A less complex autonomous soaring autopilot developed by Staffan Kjerrstrom (ref. 12) has 
been used by autonomous gliders to detect and use thermals to gain altitude. The autopilot has 
been hosted on a number of model glider airframes with mixed results. During the summer of 
2005, Alan Cocconi demonstrated high-endurance flight of a small UAV by flying his 4.27-m 
(14-ft) span solar powered motor-glider continuously for 48 hours (refs. 13-15). Autonomous 
soaring was not employed on this flight. Instead, experienced pilots guided the aircraft to soar in 
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thermals. The flight showed that a flight lasting multiple days with a small electric-powered UAV 
is possible using solar power and thermal soaring. 

A research project conducted by the NASA Dryden Flight Research Center (Edwards, 
California) developed and tested autonomous soaring algorithms using the small UAV shown 
in figure 1. The single objective of the project was to flight test a set of algorithms designed to 
autonomously detect and utilize thermals for energy gain. The autonomous soaring algorithms were 
hosted on a Piccolo Plus autopilot (CloudCap Technology, Flood River, Oregon) as an outer-loop 
guidance module computed 20 times per second. An SBXC (RnR Products, Milpitas, California) 
radio-controlled cross-country glider was chosen for flight tests because of its considerable fuselage 
capacity and excellent gliding performance. The aircraft was modified to include a flight termination 
system, a Piccolo Plus autopilot, a pitot-static probe, antennas, an electric motor, and a folding 
propeller. The resulting UAV, dubbed the Cloud Swift, has a lift-to-drag ratio of approximately 28 
and a minimum sink rate of approximately 0.4 m/s as determined from flight. 



EC05-0198-08 


Figure 1. The Cloud Swift UAV. The Cloud Swift is a modified SBXC glider (RnR Products) and 
has a span of 4.27 m (14 ft) and a weight of 6.8 kg (15 lb). 
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AUTONOMOUS SOARING GUIDANCE AND CONTROU 


The top-level autonomous soaring guidance and control structure is shown in figure 2. The 
guidance and control for soaring are implemented as an outer-loop controller to the Piccolo Plus 
autopilot. Inputs to the soaring guidance and control consist of static pressure, airspeed, throttle 
command, and the latitude and longitude of the aircraft. Estimates of the total energy rate and total 
energy acceleration of the aircraft are used by the thermal identification equations to determine 
the radius, vertical velocity, and position of the thermal. The circle guidance calculates the steady- 
state turn rate, position error, and velocity error for tracking a circular path inside the thermal. The 
controller attempts to drive position and velocity errors to zero and respond to changes in energy 
acceleration. Finally, mode logic uses aircraft energy terms to determine when the aircraft should 
be searching for thermals (by flying waypoints) and when it should circle within a thermal. 



Turn rate 
command to 
Piccolo Plus 
autopilot 
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Figure 2. Top level block diagram of autonomous soaring guidance and control. 


Total Energy Estimation 


Static pressure sensor input, estimated true airspeed, and the aircraft throttle command were 
used to determine the rate and acceleration of aircraft total energy. The energy rate and acceleration 
define the basic parameters that are used to guide the aircraft during soaring flight. Aircraft total 
specific energy was first calculated by collecting the potential and kinetic energy terms of the 
aircraft energy equation normalized by aircraft weight. The resulting energy term, E , is referred to 
as total energy in this paper and was calculated using eq. (1). 


E = h 


est 



( 1 ) 
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Equation (1) is solved using gravity, g, airspeed, V, and the estimated aircraft altitude,/?^, from 
a representation of the standard atmosphere (ref. 16). Equation (2) gives a reasonable fit to the 
standard atmosphere for altitudes below 10 000 m. 


est 


f y0.1903 

P static 

101325 


44333.7 


( 2 ) 


An approximate pressure altitude estimate was adequate for autonomous soaring because only 
the rate of change of energy is used by the guidance and control. The total energy calculations are 
illustrated in figure 3. Total energy was filtered and differentiated with filter 1 and filter 2 to obtain 
specific energy rate, E . 



command Engine . 

effect on E 


Figure 3. Total energy estimation. 
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Energy rate was corrected to account for the energy added to the aircraft by the motor. The 
correction term was calculated by differentiating eq. (1) and solving for motor terms only. The 
resulting relationship is given in eq. (3). 


W 

_ motor 
^ motor — 

O 


(3) 


Equation (3) was solved using an approximation of V mot0f . given in eq. (4) 


V, 


motor 


T 
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(4) 
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where T is the axial force from the motor and m is the aircraft mass. Thrust is calculated from the 
simple engine model given in eq. (5) 


T = 


0.64(//?ro/ / _ | - throi [nm ) T n 
s 2 + 0.16^ + 0.64 


( 5 ) 


where throt. 


I is the 1-s delayed throttle command and throt trjm 


is the trim throttle needed for level 


flight. The term T is the maximum thrust generated by the motor. One second was chosen for 
the throttle delay to account for delay in the motor and speed controller systems and to approximate 
the delay in the E calculation. Equation (5) is essentially a gain of T max on the change in delayed 
throttle command delta from trim throttle and a second-order filter with a 0.8 rad/s cutoff frequency. 
Values for T and throt. 
limited at ±10 m/s 2 to reduce transient behavior. 


. max . were taken from flight data. The estimated energy rate was rate- 


Filter 3, shown in figure 3 was used to calculate energy acceleration, E . Finally, energy 
acceleration was modified to produce scaled energy acceleration, E s , using eq. (6). 


E= : 

max(F,1.0) 


(6) 


Scaling was used to reduce the sensitivity of the controller to variations in thermal strength. The 
use of scaled energy acceleration allowed the soaring controller to behave similarly while flying in 
weak lift with small variations in energy acceleration or flying in strong lift with large variations 
in energy acceleration. 

The three filters used in the total energy calculations have two purposes. First, the filters must 
remove noise from the aircraft energy acceleration estimate that is caused by imperfect pressure 
sensor accuracy and resolution. Second, the filters must twice differentiate total energy to produce 
energy acceleration. The undesired effect of this method is a high time delay. This estimation 
latency reduced the effectiveness of the E feedback channel in the controller and was the primary 
reason for the development of the thermal identification algorithms. Possible ways to reduce delay 
in future soaring algorithms include the use of better sensors, higher-order filters, or incorporation 
of aircraft acceleration into the energy state estimation equations. 

Thermal Identification 

An estimation of the thermal location, size, and strength was made to improve soaring 
performance. The top level block diagram for the thermal estimation is given in figure 4. The 
aircraft position, P , was calculated in Cartesian coordinates from the latitude and longitude using 
eqs. (7) and (8). 
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Figure 4. Thermal estimation top level block diagram. 


P x = ( Lat - Lat^a 
Py = iyLon - Lon^acos^Lat^ 


(7) 

( 8 ) 


Equations (7) and (8) are solved using the equatorial radius of the Earth, a , the aircraft latitude 
and longitude, Lat and Lon , and the latitude and longitude of the chosen origin location, Lat 0 
and LoUq (such as the location of the UAV ground station). 

A history of previous aircraft positions and energy rates was collected in a first-in-first-out 
queue, given by eq. (9). 
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(9) 


Entries to the queue given in eq. (9) were made using a queue sample time step, dt , of 1 s. 
The subscript k refers to the measurement index. The length of the queue, M , used in eq. (9) was 
chosen to be 45 samples because the aircraft roughly completes 2 circles in tliat time. 
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Estimation of the movement or drift of the thermal was found through flight test to be 
important for achieving a good estimate of thermal position and radius. Drift is primarily caused 
by wind, but early flight tests showed that thermals do not always drift with the same velocity or 
direction as the prevailing winds. A discussion on thermal drift with possible reasons why it is not 
equal to the wind velocity is given in reference 17. Thermal drift was calculated by comparing the 
position of the thermal given by the first M D entries of the queue with the position of the thermal 
given by the last M D entries of the queue. 

First, the energy rate entries of the queue were shifted using eqs. (10a) and (10b) to prevent 
negative values from entering subsequent calculations. 

^(vE 1,2 ]) = ^i(v[ 1,2 ]) (10a) 

^(v 3 ) = ^(v 3 ) _min (^( : ’ 3 )) (10b) 

Equation (10) produces the shifted queue, q s , where the term q^ (:,3) represents the third column 
of <7| . The thermal position given by the first or most recent M D rows of q s was calculated using 
eq. (11). 
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The thermal position given by the last M D rows of q s was calculated using eq. (12). 
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( 12 ) 


The estimated drift velocity of the thermal, V D , was calculated from eq. (13) 

P, - R, 


V D = 


thD 1 thD 2 


M - M n * dt 

q D) q 


(13) 
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The resulting drift velocity was limited to ±10 m/s and rate limited to ±0.1 m/s 2 . A restrictive rate 
limit was chosen to retain previous drift calculation information. The estimated drift velocity was 
then used to form a drift-corrected queue, q 2 , using eqs. (14a), (14b), and (14c). 


= (Z ? " V V d( 1 ) * dt q + ^V 1 ) 

(14a) 

<?2 (V 2) = ~ X ^ V d{ 2 ) * dt q + ?j(V 2 ) 

(14b) 

^(v 3)= ^ ( v 3) 

(14c) 


Figure 5 shows the mapping of aircraft position values in q s to the positions given in q 2 to 
account for drift. Typical soaring flight data was used in this illustration. The result of the drift 
correction is to transform the aircraft position measurements into a reference frame that is moving 
with the thermal. The resulting queue, q ~, , can then be used to give a better real-time estimate of 
the thermal position and radius. 
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Figure 5. Comparison of the queue q s , before the drift correction, to the queue q c , after the drift 
correction. Drift was calculated using P t j lD ^ and P r j lDi ■ 
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The estimated position of the thermal center, P ., , was calculated using the corrected queue, 
q 2 , as given in eq. (15). 


M q M q 
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(15) 


Equation (15) represents the position centroid of the aircraft energy rate during the previous M 
samples. This method for estimating thermal position is easy to implement, does not require 
tuning, can be used for any thermal size, and does not have a high computational cost. Drawbacks 
associated with this method are sensitivity to the length of the queue and a bias of the estimated 
thermal position toward the mean position of the measurement set. The centroid calculations found 
in eq. (15) differ from the centroid calculations found in eqs. (11) and (12) because the energy rate 
term in the numerator is squared in eq. (15) to increase the sensitivity of the estimated thermal 
position to changes in energy rate. 


Thermal vertical velocity, w fh , was approximated from the aircraft energy rate history found 
in g | using eq. (16). 


w th = l.lmax^O, 3 )) (16) 

Equation (16) uses the assumption that the maximum thermal vertical velocity is 10 percent greater 
than the maximum energy rate experienced by the vehicle during the previous M samples. This 
very simple estimation method was used because high accuracy was not required. The thermal 
vertical velocity was rate limited to +0.025 m/s 2 and -0.015 m/s 2 as a method to retain old 
information about the thermal strength that has been passed from the queue. The estimated thermal 
vertical velocity was used in addition to the values in q 9 to estimate the thermal radius. 

Thermal Radius Estimation 

The radius of the thermal was estimated using a gradient descent method to fit an assumed 
thermal shape to the data found in the queue. The estimation process started with the calculation 
of the distance between each aircraft position in q 2 to the estimated thermal center, P th , using 
eq. (17). 


S(i q ) = “ P thO)) + (^ 2 (^ 2 ) - P th ( 2)f 


(17) 
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The vector, S , was used to calculate the predicted thermal vertical velocity of each position in 
q 2 using the thermal vertical velocity distribution developed by Wharington (ref. 11) but modified 
here to account for the environment sink rate, V e . The resulting shape is given in eq. (18). 


W pred 





max(r t h , 1 ) 

V J 


V„ 


(18) 


Equation (18) was first solved using an initial value of 45 m for the estimated thermal radius, 
r f j : , and estimated thermal radius from the previous frame thereafter. The shape of the predicted 
thermal vertical velocity distribution is plotted in figure 6. The radius of this shape is located at 
w pred / w th ~ 0-24, not at the zero crossing. 
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Figure 6. Predicted thermal vertical velocity distribution from eq. (19) for w fh = 2 m/s and 


w, 


th 


45 m. 


A weighted error, e w , was calculated by taking the difference between the predicted thermal 
vertical velocity, w ^ , and the energy rate of the aircraft given in the third column of q 0 . The 
weighted velocity error was calculated using eq. (19). 

e J l q) = [ W pred^q) " (*<,>*)) H q 0 q ) ( 19 ) 

The velocity influence vector, H , was used to weight the more recent values of <? 2 and w prec j 
higher than the older values. Values for H are plotted in figure 7. 
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Figure 7. Influence vector, H , for a queue size, M , of 45. The influence vector was used to 
de-emphasize older data points. 


The mean sum squared velocity error, e , was calculated using eq. (20). 

o liiHo U 
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sumsq 
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( 20 ) 


Equation (20) gives a measure of the fit between the assumed thermal shape given in eq. (18) 
and the measured energy rate of the aircraft stored in the queue. The sum squared velocity error for 
a perturbed radius, e sumsq , was calculated in a similar way using eqs. (21) -(23). 



( 21 ) 

( 22 ) 


e sumsq 


e w * e w 




(23) 
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The perturbed radius, r., , used in eq. (21) is defined by eq. (24). 


' th 


= r th +5 r 


( 24 ) 


The stepsize, <5 ; . , was chosen to be 0.5 m for this study because this value gave smooth and rapid 
convergence of the estimated thermal radius and was robust to thermal size variations. 


The error for the current thermal radius and the error for the perturbed thermal radius were 
used to calculate a gradient as given in eq. (25). 
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p — p 

sumsq sumsq 


(25) 


The gradient was then used to calculate the thermal radius for the next frame using eq. (26) 


'th 


' k+\ 


= L.A.. + r, 


th 


(26) 


with a learning rate, L y , of 10. A learning rate of 10 was chosen because it gave rapid and smooth 
convergence of the estimated thermal radius with flight data. Estimated thermal radius was limited 
to be greater than 40 m to prevent soaring in small diameter thermals requiring high bank angle 
flight and to be less than 80 m to make the soaring algorithms more robust to variations of the 
position and strength of large thermals. Alternative methods for estimating thermal radius exist 
(ref. 11) but were not used here. 


Guidance Calculations 


Guidance calculations were used to determine the position error, velocity error, and the 
steady-state turn rate associated with flying in a circular path inside the thermal. The aircraft was 
commanded to fly a circular flight path having a radius, r ^ , determined by eq. (27). 


r cmd = °- 65r th 


(27) 


The scale factor 0.65 was chosen because it defines a flight path radius that is usually small enough 
to be within the thermal “core” but large enough to avoid the high aircraft sink rate associated with 
high bank angle turns. The steady-state turn rate, \j/ , needed to fly at an airspeed of V along a 
circular path given by r cmd is given in eq. (28). 


Vss = 
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V cmd / 


(28) 
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The negative sign in eq. (28) ensured that the aircraft always circled to the left while soaring. 
This was done to give the research personnel a quick visual indication of the aircraft mode during 
flight test. Flight plans were designed to command right turns when the aircraft was searching for 
thermals and left turns when the aircraft was in soaring mode. Future soaring guidance algorithms 
could potentially reduce the time needed to initially center a thermal by using the aircraft roll 
response to determine the turn direction as described in reference 1. 

Controller position error, e p , was generated using eq. (29). 

= r cmd ~ SO) ( 29 ) 

Position error was differentiated and filtered to create controller velocity error, e , using the transfer 
function given in eq. (30). 


e 


V 


1.695 

5 2 + 2.65 +1.69 


(30) 


The guidance calculations given in eqs. (27) -(30) calculate the position error, velocity error, and 
steady-state turn rate information that was used by the soaring controller to maintain a circular 
flight path centered on the thermal. 


Soaring Controller 

The soaring controller architecture is loosely based on one of the thermal centering methods 
given for glider pilots by Reichmann (ref. 7). A paraphrase of Reichmann’s method gives three 
rules: 

1. As climb improves, flatten the circle (approximately 15-20° bank angle) 

2. As climb deteriorates, steepen the circle (approximately 50° bank angle) 

3. If climb remains constant, keep constant bank (approximately 25-30° bank angle) 

These rules were used as a starting point for the soaring controller design shown in figure 8. Energy 
acceleration feedback was used to implement rules 1 and 2 because the rules given by Reichmann 
are to be applied to the “changes in rate of climb.” The steady-state turn rate term kept the aircraft 
at a constant bank angle when the feedback terms were zero to fulfill mle 3. The soaring controller 
inputs consist of scaled energy acceleration, E s , position error, e p , velocity error, e v , and steady 
state turn rate, V f ss . The soaring controller output consists of a turn rate command, V C md , that 
was sent to the autopilot. The energy acceleration gain, , was tuned in simulation to produce 
an aircraft response that followed the rules outlined by Reichmann. The thermal model given in 
reference 18 was implemented into a hardware-in-the-loop simulation for controller evaluation 
and gain tuning. The k j and k 2 gains were tuned to make the path of the aircraft more circular 
and to improve the damping of the controller during soaring flight in a variety of thermal sizes 
and strengths. The resulting control gains were chosen to be 50, 0.4, and 0.165 for k ^ , k 2 , and k 2 
respectively. 
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Figure 8. Soaring controller. Turn rate command output from soaring controller is tracked by the 
autopilot during soaring flight. 


Mode Logic 

Engagement and disengagement of the soaring controller during flight was handled 
autonomously using the mode logic shown in figure 9. The mode logic described in this section 
is significant because it determines how a thermal is detected and what conditions must exist 
before the aircraft will leave the thermal. The mode logic is broken into three components; thermal 
detection, thermal rejection, and a latch mechanism. The thermal detection logic uses aircraft 
response to detect thermals. An encounter with a thermal by the aircraft produces large changes in 
energy rate and energy acceleration. As the airplane flies a straight path through a typical thermal, 
the energy acceleration reverses sign and becomes negative when the aircraft is nearest to the 
center of the thermal. The thermal detection logic is tuned to trigger when the thermal is imputing 
significant energy to the aircraft and the aircraft is closest to the thermal center. 

The decision to switch from soaring flight to searching flight was handled by thermal rejection 
logic using energy rate. Thermal rejection was needed to switch from soaring flight to searching 
flight when the thermal was no longer providing sufficient energy to the aircraft. Two energy-rate 
based criteria are used by the thermal rejection logic. First, thermals are rejected if the output of 
filter 4 is less than -0.5 m/s. This criterion would be met if the aircraft encountered light sink for 
a period of approximately 5 s. The second criterion, based on the output of filter 5, triggered if the 
aircraft had not experienced sufficient energy rate for approximately 10 s. Additional criteria, not 
shown in figure 9, were used to disengage from soaring mode when the upper or lower altitude 
limit was exceeded, a loss of communication between the UAV and ground station (lost link) 
occurred, or soaring mode was disabled by the ground station operator. 
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Figure 9. Mode logic diagram. The aircraft will enter soaring flight when Soarmode = 1. 


RESULTS 

A total of 17 flight tests with the Cloud Swift were conducted during the summer of 2005 
to perform aircraft checkout, tune autopilot gains, and test the autonomous soaring guidance and 
control. Twenty-three thermals were detected and utilized during the soaring research flights. The 
average altitude gain obtained from a single thermal was 173 m (567 ft) as measured by the Piccolo 
Plus GPS. 

The path of the aircraft during a typical soaring flight in an average thermal is shown in 
figure 10. The aircraft encounters two thermals during this time period. A small thermal at position 
P x = -900 and P y = 300, was not sufficient to cause the mode logic to engage soaring mode. The 
second thermal encountered at P x = -1250 and P y = 100 barely imparted sufficient energy to the 
aircraft to cause a mode change. The estimated thermal position, shown in black, moves as the 
aircraft encounters increasing or decreasing energy from the thermal. During soaring flight, the 
aircraft progresses in the positive X and Y direction (North and East) as the thermal drifts. The 
aircraft attempts to increase the time it spends in lifting air as it circles in a thermal that is likely 
varying with time and altitude. 

Figure 11 shows the energy state of the aircraft for the flight path shown in figure 10. 
Controller performance was determined from the time history of energy rate. Energy rate should 
be constant when the aircraft is perfectly centered in the thermal. In this time history, energy 
rate fluctuates as the aircraft climbs in the thermal, indicating that improvements in the soaring 
algorithm can be obtained. Postflight estimated energy rate E j - , shown in the second subplot of 
figure 1 1, represents the true energy rate of the aircraft and was obtained by first forward-backward 
filtering the aircraft energy with a second-order low-pass filter and then differentiating the result. 
The onboard estimate of E differs from the postflight estimated E because the onboard estimate 
includes the delays caused by filter 1 and filter 2. The delay of the onboard estimate of E and E 
resulted in reduced controller performance. 
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Figure 10. Aircraft path during soaring flight in an average thermal. The aircraft starts at the upper 
right and flies in a straight line until a sufficiently strong thermal is detected at the bottom left 
comer of the figure. 
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Figure 11. Aircraft energy state during soaring flight in an average thermal. 



Figure 1 1 also shows the throttle command drop to zero shortly after the thermal is encountered, 
indicating gliding flight with the propeller folded. The Piccolo Plus autopilot was configured to use 
the throttle to control altitude. This caused the throttle to go to zero and the propeller to fold back 
as the aircraft climbed above the commanded altitude while soaring. 

The bank angle, pitch rate, and yaw rate are shown in figure 12 for the section of flight shown 
previously in figures 10 and 11. Variations in bank angle reveal the activity of the aircraft as it 
attempts to stay within the thermal. 



Soaring mode 

switched on 


Figure 12. Aircraft states during soaring flight in an average thermal. 
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Thermal position estimation values are shown in figure 13 for the same thermal. Aircraft 
path and drift corrected path were taken from and g 9 respectively at 232 s. The figure shows 
increasing thermal velocity toward the right side of the drift corrected path causing the estimated 
thermal center to move in the positive direction. The distance of each position in q 1 to the estimated 
thermal center is plotted against the energy rate of each point in figure 14. The estimated thermal 
shape is also plotted in figure 14 to illustrate the performance of the thermal estimation algorithms. 
The figure shows that the basic shape of a thermal can be approximated by eq. (18); however, the 
estimated thermal shape did not fit the data perfectly, indicating that eq. (18) did not capture all of 
the properties of the thermal. 
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Figure 13. Thermal position estimation values for soaring flight in an average thermal. Aircraft 
path and drift corrected path values were taken from q j and q 2 respectively. 
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Figure 14. Thermal radius estimation values for flight in an average thermal. Estimated thermal 
velocity and radius is 2.52 m/s and 60 m respectively. 
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Figure 15 shows a failure of the UAV to detect and circle within a weak thermal. The mode 
logic in this case did not switch to soaring mode quick enough to allow the aircraft to stay within 
the lift. The switching delay was caused by two factors: the non- symmetrical shape of the thermal 
vertical velocity as the aircraft flew through the thermal and the delay in the energy rate and energy 
acceleration estimates. This problem can likely be remedied through better tuning of the mode 
logic and soaring controller and the use of an energy-state solution with less delay. 
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Figure 15. Example of a late mode switch into soaring mode. 
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The soaring controller was able to keep the energy rate more constant when soaring in a 
stronger thermal as shown in figure 16. The aircraft gained over 600 m (2000 ft) in this thermal, as 
reported by GPS, and experienced an average climb rate of 2.8 m/s (5.4 kn). The estimated thermal 
radius was saturated on the predefined upper limit of 80 m during this climb. Energy rate and 
acceleration delay may have had less impact on controller performance as the aircraft flew larger 
circles in stronger lift. Figure 17 shows the altitude gain along with a more constant E associated 
with improved controller performance. 
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Figure 16. Aircraft path and estimated thermal center during soaring flight in strong lift. 
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Figure 17. Aircraft energy state during soaring flight in strong lift. 
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Although long duration flight was not an objective of this project, a partial demonstration 
of increased endurance is shown in figure 18. The top plot in figure 18 shows the altitude of the 
aircraft during one of the research flights. The cruise altitude for this flight was 960 m. The bottom 
plot shows throttle command where a value of 1 is full throttle, used for takeoff, and a value of 0 
is power off, gliding flight. The altitude time history of the aircraft during this flight is similar to 
those predicted by simulation (ref. 6) and to those of migratory birds (ref. 4). The flight was ended 
when the capacity of the actuator batteries was consumed. Despite this limitation, the UAV flew 
over 60 min without using the motor during this flight. 




Time, min 
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Figure 18. Extended soaring flight. Multiple thermals were used to glide under autonomous control 
for over 60 min. The flight was ended prematurely because the vehicle was not configured with 
actuator batteries for long endurance flight. 


CONCLUDING REMARKS 

A guidance and control method for autonomous soaring flight in thermal updrafts was 
developed and tested using a small electric-powered unmanned aerial vehicle. The total energy rate 
of the vehicle was estimated from static and total pressure measurements using digital filters and 
differentiators. The position, radius, and strength of the thermals were estimated using a history of 
total energy rate and aircraft position measurements contained in a first-in-first-out queue. Thermal 
drift velocity was estimated and used to map the aircraft position data in the queue to a reference 
frame attached to the thermal. The data in the queue was then used to estimate the thermal center 
location using a centroid calculation. The radius of the thermal was estimated using an iterative fit 
to an assumed thermal vertical velocity profile. 


23 



Aircraft velocity command, thermal position, and thermal radius were used to generate a 
circular flight path command for soaring. Position and velocity error to the flight path command 
were used along with calculated steady- state turn rate to generate a turn rate command in the 
soaring controller. The general structure and initial gain values of the soaring controller were taken 
from piloted soaring rules. Lastly, simple mode logic based on the aircraft energy state was used 
to engage and disengage the soaring controller. 

The soaring algorithms were flight tested at Edwards Air Force Base using a 14-ft span 
electric -powered UAV. Results show that the UAV could autonomously detect and exploit thermals 
to gain altitude and conserve energy. Controller performance was determined by inspecting the 
unsteadiness of the aircraft energy rate during soaring flight. Controller performance was found 
to be better in strong thermals. Latency in energy rate and energy acceleration estimates was 
found to be the primary cause of reduced controller performance in weak thermals. Future soaring 
algorithms should use better pressure sensors or incorporate aircraft acceleration information to 
improve the energy state estimate. Although demonstration of long endurance flight was not an 
objective of this project, multiple thermals were used in a single flight to soar autonomously for 
over 60 min. Overall, the flight results show successful autonomous soaring of a small UAV using 
the algorithms given in this report. 
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